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Abstract 

We derive a general formalism for constructing a four- dimensional effective action for 
the radion field in 5 dimensions, taking into account the possible dependence of the scale 
factor on the internal coordinate. First and second variations of the action with respect to 
the size of the extra dimension reveals two types of algebraic constraints on the combination 
of brane and bulk Lagrangians that extremize the action, and provide an effective mass for 
the radion field. We derive a compact formula for the radion mass squared in terms of 
brane and bulk potentials for the case of a bulk scalar field. We further specialize these 
constraints for several, five- dimensional realizations of brane world models. 



1 Introduction 



Theoretical models with internal compactified space-like dimensions play a special role 
in string theory 0, particle physics @, and cosmology @, |J ||. To be consistent with 
phenomenological requirements, the volumes of internal dimensions have to be stabilized, 
i.e. the corresponding moduli, or radions, have to acquire a regular mass term. Depending 
on the model at hand, a mass for a radion field may vary in a wild range from 10 -3 eV 
to the Planck scale. The stabilization of extra dimensions is a phenomenological necessity 
both in models with a smooth distribution of matter in the extra dimensions and in models 
with branes, where the metric may acquire a strong dependence on the coordinates of 
the extra dimensions, a = a(y,z, ..). One of the most interesting recent developments for 
model-building is the idea that the strong dependence of the scale factor on the coordinates 
in extra dimensions may have profound consequences for the gauge hierarchy problem in 
particle physics [|], [|. 

The presence of matter in these models, either placed on the brane or smoothly dis- 
tributed in the bulk, usually has a linear coupling to the radion field, displacing the radion 
from its vacuum value, and creating the T™ pressure-like components of the energy-momen- 
tum tensor, where n labels the internal dimensions. The presence of this extra-dimensional 
pressure enables one to recover the usual 4D cosmological expansion [^] and cures the crisis 
in "brane cosmology" that consists of a drastically different cosmological evolution at 
late times. The firm connection between the stabilization mechanism and extra-dimensional 
pressure was further exemplified in Refs. || |§. 

The static solution for the radion field equations (b = 0) implies that the energy density, 
"normal" 4D pressure of matter and the extra-dimensional pressure must be related 



by an integral condition. The first of a series of conditions was found in Ref. [pTQ| , in the 
context of the Randall-Sundrum model, and presented as a topological constraint necessary 
for the consistency of the particular solution. A similar integral condition, suitable for a 
non-warped case, a(y) = const, was found in Ref. and it was shown to be associated 
with the stabilization of the size of extra dimension. Finally, a whole family of different 
sum rules was presented in Ref. |H||, applicable to the case of matter that is uniformly 



distributed (^^-independent) in "our" dimensions. One of these sum rules was shown [12 



to provide a consistency check for solutions with a constant radion field (b = 0). The same 



sum rules were generalized in 6 dimensions in |13| and in d dimensions in ||14|| . 

This paper extends the analysis given in Ref. |§ and includes the case where the scale 
factor depend on the extra-dimensional coordinate of a 5D spacetime. This results in a 
constraint on the components of the higher-dimensional energy momentum tensor, the ex- 
tremization constraint, that generalizes the simplest constraint obtained in |S]]. Using the 
fact that the variation over b with branes at fixed positions, y\ and y2, is equivalent in 
fact to the variation over the length of the extra-dimensional circle, we reduce integral con- 
straints to a set of algebraic relations between brane Lagrangians and the bulk Lagrangian, 
calculated at the position of the branes. 

The presence of a stabilization mechanism leads to the radion relaxation to its minimum 
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at b = bo, with b = 0. However, the inverse is not true. The existence of a solution with 
6 = might be a result of a fine-tuning that does not lead to a physically viable "brane 
model" . Indeed, small perturbations of the solution around the extremal point b = bo may 
prove to be unstable. Such an instability can manifest itself as a tachyonic mass and/or 
a ghost-like signature for the kinetic term of the radion. To decide whether a particular 
brane world scenario is stable against radial perturbations or not, one has to find a close 
family of solutions around b = bo and expand in b — bo . In many existing brane- world 
scenarios such an analysis is possible, whereas in some models that specify a(y) and 
only at a given point b = b , the stability cannot be analyzed. In cases when it is possible, 
the analysis of the stability (b/b < 0) can be done with the use of Einstein's equations 
T5| or directly in the action |16|, [IT]]. We choose the latter method, and obtain a rather 



compact algebraic formula for the radion mass squared, again in terms of brane parameters 
and bulk quantities calculated at the position of the branes. The stabilization constraint 
can be regarded as the condition on the positivity of the radion mass squared. 

After obtaining the generic forms of the 6 = and b/b < constraints, we apply them 
for some particular realizations of brane-world models. We exemplify how the extremization 
constraint (b = 0) is satisfied in a number of static, five-dimensional models, namely, in the 
Randall- Sundrum model in the two-brane model with a Casimir force |18|, [l9| , and in 
the two-brane model with a classical massless bulk scalar [20fl. In the above cases, we also 
analyze the stabilization constraint and determine whether, or under which assumptions, 
the radion has a positive or negative mass squared. The results of this paper can be used 
to check the consistency of various brane model solutions and present a useful formalism 
that helps to analyze the stability of a particular extra-dimensional solution against small 
perturbations of the radius of the internal spacetime. 

This paper is organized as follows. Section 2 presents the complete study of the 5D 
case addressing the generic form of the stabilization and extremization constraints. Section 
3 applies the derived constraints to particular brane-world models. Section 4 presents our 
conclusions. 



2 The radion potential in 5-dimensions 

Here, we assume that the brane self-energies, Vi, do not depend on the four-dimensional 
coordinates either explicitly or implicitly. Therefore, V^'s can be, at most, functions of 
the extra coordinate y. In this case, we may write down the following ansatz for the 
five-dimensional spacetime 

ds 2 = a(t, yf{-dt 2 + c% dx i dx j ) + b{tfdy 2 (2.1) 

where a(t,y) stands for the conformal factor multiplying a four- dimensional, flat line ele- 
ment and b(t) denotes the scale factor of the extra dimension. 
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We first focus on the gravitational part of the theory. For the above ansatz, the 5D 
scalar curvature takes the form 

n = y) - + «« = | - + (-^ - %) . (2.2, 

The last two terms inside brackets, denoted as RS y \ are the part of the five- dimensional 
scalar curvature involving solely derivatives with respect to the extra coordinate y. The 
covariant derivative of the scale factor b(t) will turn out to be a total derivative in 4D 
spacetime, and hence drops out. The scalar curvature R^ft^y) is still ^/-dependent and, 
therefore, the purely time-dependent, four- dimensional scalar curvature needs to be 
extracted. For this reason, we define a(t,y = 0) = ao(t) and write the conformal factor as 

a(t,y) = a (t)a(t,y) (2.3) 

where a(t, y = 0) = 1, by definition?]. a (t) will play the role of the 4D metric component 
that we need for the effective theory and so we write 

R^(t, y) = 1 ^ (4) - 4t d»a] , (2.4) 



or a°a. 



o 



with = 6ao/ag 



i . 

We now turn to the action of the five- dimensional gravitational theory, which can be 
written as 

R 



S 



J d 4 x<fdyV^G{-^ + c} (2.5) 

where k\ stands for the five- dimensional Newton's constant and £ is a generic Lagrangian 
that may include bulk cosmological constants, brane-self energies and various bulk and 
brane matter fields. The integral over y is performed over a compact dimension. We can 
recover finite volume, but infinite extra dimension type models by sending one of the branes 
off to infinity. In case when two infinitely thin branes are located at y% and y2, the domain of 
integration over y can be decomposed into three separate pieces that represent the integral 
over the bulk and the integrals across the branes: 

V2 ryi ryi+£ ry2+e 

+ / +/ +/ • (2-6) 

yi Jy2 J y\-t Jy2—t 

The part of the above action consisting of the five- dimensional scalar curvature, R, will 
provide the kinetic part of the four-dimensional effective action describing the dynamics of 
the graviton and the radion field. By using the expressions ( p.2[) and (|2.4j) and integrating 
by parts, the gravitational part takes the form 



Sg = J d A x <j> dy ^/—G 



R 

24 



J rf 4 xv^| ^dyd 2 b^^ + ^ j dyd li {ab)d»a + f dyd 4 b^^, (2.7) 



1 Note that, in the case of static branes, this definition is not necessary as ao(t) reduces to a constant. 
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with \J — G = <2q a 4 b and \f--gl 



In order to obtain standard 4D Einstein gravity 



we perform a conformal transformation of the four- dimensional metric tensor in order to 
eliminate the coupling between the radion field and the scalar curvature TZ^ A > . We set: 



(9 



4)fiu — 

where 

which leads to the result 

K {A) = A(b) 



A{b) 



is- 



A 2 (b) 



A(b) 



n 2 b 



dy a (t, y) , 



n {4) + ZD^D* In A(b) <9 M In Ad^hxA 

2 



(2.1 



(2.9) 



(2.10) 



Then, the entire action can be written as, 



S = - 



J d A Xy/^gl 



o2 +B(b)d,bd»b + V eff (b) 



with 



and 



(/ dy d^ab) d^a) 3 ,9 In An 2 



d^bd^b 



db 



Veff(b) 



A 2 (b) 



dy a b 



R(y) 
2^1" 



(2.1i; 



(2.12) 



(2.13) 



Note that, having restored the four- dimensional gravitational term, all remaining terms in 
the action will constitute the building blocks of the effective action for the radion field. 
Therefore, terms involving x^-derivatives of the auxiliary functions a and A(b) are bound 
to contribute to the kinetic term of the only remaining degree of freedom in the theory, 
that is of the radion field, while any other terms will contribute to its effective potential. 

At this point, we are able to obtain the constraints on the combination of the scale 
factors and components of the stress-energy that were previously derived in || and [|Il[ . 
These constraints were derived under the assumption that 5d/5b = 0, that is, that the 
perturbation with respect to the value of the radion field does not affect the value of the 
scale factor. In this case, A(b) is simply proportional to b and the variation of the action 
3) takes the following simple form: 



5S 



d xdy 
d 4 x dy 



5(V^GR) 6G 



AIN 



oGmn 
2k\ bA 2 
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•ia 
a 



+ 



2 

12a' 2 



1 oG 



MN 



a- 



3.14) 



— <j> d A x dy \/—G 



12a" 
K 2 b 2 a 



+ T£ - 2T| 



5b 
2b 
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Note that the variation in the quantity b 2 R vanishes with these assumptions. In the above 
expression we performed an integration by parts and used the standard definition for the 
stress energy tensor. It is easy to see that the square brackets in the second line of Eq. 
(|2.15| ) simply contains a combination of the ii and 55 components of Einstein's equations, 
Rff — ^G^Rx — K^T^f = 0. Obviously, Eq. fl2.15|) must vanish identically when a(y) is the 



solution. Integrating by parts the term with a" we obtain the following constraint: 

"36a /2 



<j> d A x dy \f—G 



(2.15) 



This constraint generalizes a formula from and includes the dependence of a on y (i.e. 
warping). The constraint derived in |J is valid only in the case of small Tff. Indeed, a non- 
trivial dependence on y-coordinate usually arises due to the presence of branes. In this case, 
a 1 /a is proportional to the brane energy density p. In this case, the first term in ( |2.15| ) is just 



an 0(p 2 ) correction and can be dropped, leading exactly to the expression found in ||. For 
5-dimensional spacetimes with significant warping, the constraint ( |2.15[ ) is valid instead. 



A similar, but not identical, constraint was derived in Ref. [|yj for solutions with large 
warping by rearranging the components of Einstein's equations. The same combination 
— 2T| appears in both constraints, however Eq. ( [2.1 5|) has an extra term, proportional 



to a' 2 /a 2 , and a different coefficient in front coming from \J—G. The apparent differences 
are due to the fact that the two constraints were derived with different methods, however, 
they are both valid under the same assumptions and, are therefore, equivalent. 

In the above derivation we assumed that a and b are totally independent, and this is 
why the form of the constraint comes out to be identical with the combination suggested 
by Einstein's equations. However, in order to study the stability of the physical size of the 
extra dimension, one can no longer assume that a and b are independent and that A(b) is 
simply linear in b. Instead, we would have to consider a as a function of by, a = d(by) and 
the variation of action over b might be a more complicated procedure. It is important to 
note that there is always a residual symmetry under a simultaneous rescaling of y and b. 
All physical parameters must remain invariant under y — > Cy and b — > b/C, where C is an 
arbitrary constant. Below, we will consider a two-brane model with branes positioned at 
yi and y 2 - Then, the extremal size of the extra dimension is given by certain fixed values of 
yibo and 3/2^0, which are invariant under the above symmetry. The variation over b around 
bo with fixed y\ and y 2 is equivalent to a variation over the length of the fifth dimension. 
Moreover, using this rescaling symmetry, we can always choose bo — 1; i- e - absorb b into 
the definition of y\ and y 2 - 

A canonically normalized radion field r(t) is obtained by demanding that the kinetic 
term 2B(b) d^bd^b can be written as d^rd^r. For small deviations of the scale factor b(t) 
from a stable minimum, bo, of its potential, the normalized radion field may be written as: 
r{t) = ^2B(b ) [b(t) — bo), if the radion is not a ghost-like degree of freedom, i.e. B > 0. 
Note that, in the pursuit of deriving extremization and stabilization constraints, we can 
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continue working with the metric function b(t) instead of the radion field r(t) since: 



dV e 



eff 



dr 



,db^ dV f 



eff 



\Qr' db 



d 2 v f 



eff 



bo 



dr 2 



,dbx2 d 2 V c 



eff 



ro 



^dr' db 2 



bo 



(2.16) 



Therefore, for a well defined relation between r(t) and b(t), the transition from one field to 
the other makes no difference in the vanishing of the first derivative of the effective potential. 
In the same way, stability of the extra dimension in terms of b(t), that is d 2 V e ff > 0, leads 
to stability in terms of r(t), d 2 V e ff > 0, and vice versa. We note that in particular models, 
the value of B(bo) may also be very important because it leads to a rescaling of the mass 
term, and a very large value of B(bo) may lead to a phenomenologically unacceptably low 
value of the radion mass. 

Before we go into explicit calculation of the variation of the action over b with in- 
dependent a, we would like to specify some components of the matter Lagrangian that 
we consider here. We allow C to include the bulk Lagrangian and brane self-energies that 
may depend on the value of bulk fields: 

/ry2 ryi 
dya 4 b£= dy~a%C B + dy ~a% C B + £ a^by^Vi + U(by u by 2 ) . (2.17) 



i=l,2 



In this expression, the bulk Lagrangian may include a variety of fields, however, for the 
sake of simplicity we will restrict to the case of a single scalar field 



(2.18) 



Note that the bulk cosmo logical constant A# has been absorbed into the definition of Vb- 
We also allow the brane self-energies to depend on the value of the same field 0. A constant 
part of Vi at the extremal value = 0(&o2/«) is a usual brane self-energy, Aj = Vi(b yi). Like 
a, <p is always a function of by, Lb = Cs(by) and Vi = Viipyi). 

Finally, U(byi,by 2 ) in fl2.17|) represents an effective interaction that may arise due to 
quantum effects in the final volume and which often cannot be formalized in the language 
of the Lagrange density in y-space. Sometimes this piece of the potential by itself may 
lead to the stabilization of the extra dimension. For the remainder of the calculations we 
would like to keep a generic form for this potential, specifying its form only in particular 
applications. 

We now turn to the explicit calculation of the variation of action with respect to b. 
In order to do that, we perform a change of variables, ^ — by. This change brings about 
a significant simplification, as b now only enters in the calculation through the limits of 
integration, i.e. at specific points y\ and y 2 



A 2 (b yi ,by 2 )V e 



eff 



rby 2 

2 / ~a\0 

J by! 



1 125^(0 
2k 2 a 2 (0 



d£+ E ~a A (by i )V i + U(by 1 ,by 2 ). 

(2.19) 



1=1,2 
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where we have now assumed equally spaced branes on the circle. In the above expression, 
and henceforth, primes denote differentiation with respect to £. Note, that we used the 
integration by part to remove a" from the R y -part of the effective potential. This allows us 
to remove the so-called Gibbons- Hawking boundary terms, i.e. singular terms in when 
a" is taken at the positions of the branes. When b is taken at its extremum, V e ff can be 
interpreted as the effective four- dimensional cosmological constant. 

The first derivative over b is computed trivially, and the result takes the following form: 



A 2 (byi,by 2 



dV f 



eff 



db 



4«2 



T Ayia (byi)V eff 

K 5 



i=2 



i=i 



12a' 2 

+ ^B(by, 



i=2 



+ E Vi 



d 



i=l i=l,2 



d{byi) 



(2.20) 



Vidf + U 



where = aQnji). The terms with VJ's can be further expanded as 







d(byi 



-Via 



~<1 

a, 



a,- a® 



(2.21; 



In this expression, the signs of a\ and d>[ are defined according to the following rule: 

da 

da 
di 



a' 



a' 



(61/1+0) 



(62/2-O) 



(2.22) 
(2.23) 



with similar definitions for d)'. In other words, the derivatives are taken on the right side 
of brane 1, and on the left side of brane 2 since we assumed y% < y%. 

To derive the extremization constraint, we put b — bo — 1, and use equations of motion 
for the scale factor and 0, as well as junction conditions on the branes, both for the scale 
factor and the scalar field: 



a' 



a 



O 

l dV 1 
~ Yd6 



~a! 



Hv2) 

a 



<t>\V2 



6 



-V 2 , 

l dV 2 
YdS 



(2.24) 



The final result takes the following form: 



1 (dVi" 



i=2 



+ E y 



d 



i=i 



1,2 d ( b Vi 



-U(by 1 ,by 2 ) 



0. (2.25) 



b=b =i 



Here we used the smallness of the four-dimensional effective cosmological constant and 
put the first term in ( 2.20|) to zero. This algebraic constraint (as opposed to an integral 
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condition (gjj) ) can be used to check the consistency of particular solutions. We note that 
in the particular cases that we study later, the separate cancellations occur for yi and 1/2 
proportional terms. 

It is easy to see that the extremization constraint is not uniquely defined. Indeed, one 
can add to ( |2.25| ) linear combinations of T^yi) — k\V 2 /Q and C — Ve^i) — (dVi/d(p) 2 /S, 
since these combinations are equal to zero due to the 55-component of Einstein's equations 
and the boundary conditions for the scale factor and the scalar field: 



~6 



V? = T B 5 (Vi) 



(2.26) 



Using the above relations, we would like to give another legitimate form of ( p. 25 ): 

i=2 



o~4 



i=l 



E Vi 

i=l,2 







d{byi 



-U(by 1 ,by 2 



0. (2.27) 



b=b =i 



In order to obtain a radion mass term and impose a stabilization constraint, we differ- 
entiate dV e ff(b)/db over b: 



A 







db 2 








K 5 



,4y 



i=2 



If 



i=l 



i=l,2 



i=1.2 



diby^dibyj 



(2.28) 



Not shown here are the terms proportional to the first and second derivatives of A(b). 
These terms exist in principle, but vanish when b = b because of A e jf = and because of 
the extremization constraint ( [2.25 ), and thus have no relevance for the calculation of the 
radion mass term. It is easy to see that the above expression contains a". This may look 
worrisome as a" has a singularity right at the position of the branes. However, in view of 
the definitions (|2.22| ) the double derivatives in Q2.28j ) should be understood as the limit of 
the continuous part of d"(y) at y — > yi. 

As in the case of the extremization constraint, the expression ( [2.28|) can be significantly 
simplified at b — bo — 1, because of the use of the boundary conditions and the equations of 
motion for the scale factor and the scalar field. The extensive use of these relations allows 
to reduce (j2.28|) to the following rather compact expression: 




E vi^ 



T^ yi ) + C B ( yi 




(2.29) 



dVi dV B 



+ E ViVj 



d 2 



i=l,2 



d{byi)d{byj 



-U 



b=i 
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The stabilization constraint is the requirement on the positivity of fl2.29| ). This result can 
be further modified to an equivalent form, using relations ( [2.26 ) and the extremization 
constraints ( 2.25 ) and ( p. 27 ). One of these forms contains only bulk and brane potentials 
and their derivatives over 6: 



A2 d 2 V eff 



db 2 

+ E» 

i=l,2 



b=b Q =i 

ld 2 Vi (dVi 



4k 2 

6 i=l,2 



(2.30) 



dVi dV B 



+ E ViVi 



d 2 



i=l,2 



d(byi)d(byj 



-U 



b=l 



The generalization of ( [2.25] ), ( |2.27| ), ( |2.29| ) and ( |2.30| ) to the case of multiple scalar fields is 
straightforward. 



3 Applications to brane-world models 



1. No branes, trivial scalar field profiles. 

In this case we take Vi = 0, = const, Vb = and o» = 1. In the absense of 
branes, the individual positions yi lose their meaning and translational invariance in y- 
direction requires that all quantities depend only on the difference 2/2 — 2/1, which is simply 
related to the size of the compact dimension. Therefore, V e ff(byi,by2) = V e ff(b(y 2 — yi))- 
Furthermore, since Vi = 0, we also have T| = (in the absence of branes this is true for all 
y). Looking at the extremization constraints, we immediately discover that 



2A B (y 2 



, dU 



0. 



(3-1) 



6=1 



The requirement A e ff = is equivalent to the condition 2A# (y 2 — y\) + U = that ensures 
the vanishing of the other components of T^. The stabilization constraint reduces to 



A 2 



,d 2 V 



eff 



db 2 



d 2 U 



b=b = l 



db 2 



> 



(3.2) 



6=1 



Particular examples of such a stabilization due to the Casimir force generated by multiple 
massive and massless scalar fields while classical bulk profiles of these fields are trivial were 



considered recently in Ref. [21 . 



2. Randall- Sundrum model: empty bulk, empty branes, U = 0. 
In this case, the extremization constraint takes the form 

i=2 



2ajy 







(3.3) 



8=1 
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and leads to the usual relation between the bulk cosmological constant and the brane tension 
in the RS model: —As = nlA 2 /6. The radion mass term, Eq. ( |2.29| ) is also zero because 
Tg + Cb — T§ + Ab — 0. It can actually be shown that every derivative of the radion 
effective potential is trivially zero for this particular solution. Going back to the expression 
of the potential, Eq. ( |2.13|) , substituting the relevant quantities inside the integral and 
performing the integration over the extra coordinate, we are led to the result 

i=2 

(3.4) 

i=i 

We therefore conclude that every derivative of the effective potential with respect to b, 
once evaluated for the specific solution, will be identically zero. Therefore, this particular 
solution is an absolute saddle point (a flat direction) of the radion effective potential. 



A 2 V, 



eff 



4\A 



B\ 



A B + 



K 



—A 2 
6 1 



3. Empty branes, trivial scalar field profile, U ^ 0. 



Next, we turn to a particular two-brane- world solution in which the Casimir force gen- 
erated by a conformally coupled scalar field was taken into account |18|, [T3|. For the case 
of an attractive Casimir force, the ^/-profile of the scale factor takes the following form: 



a(t,y) = cosh 2 / 5 (a;%) , 



UJ 



25 
24 



4\A 



B\ 



(3.5) 



For a repulsive Casimir force, the above cosh-like solution is substituted by a sinh-like 
one and the analysis follows along the same lines. The solution assumed a vanishing 4D 
cosmological constant, A e // = 0, and no bulk scalar fields contributed classically to the 
five-dimensional action. For both cases, the presence of the Casimir force is accounted for 
by the U potential in the action: 



2a 

U(byi,by 2 ) = -r— 



(3.6) 



where a is a dimensionless constant and L the inter-brane distance in 5D conformally-flat 
coordinates, defined as 



1 



I(ub) 



wby 2 



uj u Jujby! cosh 2/5 £ 

The extremization constraint ( |2.27| ) takes the following simple form: 



(3.7) 



Ab + -j-A 2 + 

6 ap 5 (uj) 



i=2 



i=i 



0. 



(3. 



Indeed, after applying the following relations from |19 

6|A fl ' 



Ai 



K 2 



tanh (uyi 



B\ 



Aau 5 



(3.9) 
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Eq. (|3.8|) is satisfied, separately for y\- and ^-proportional terms. 

Using the above relations, we can calculate the radion mass term from Eq. ( |2.30| ): 



A 2 d<2V eff 



db 2 



■10\A B \u< 



6=1 



m _ yi 

a 2 cli 



6 



v\A 2 | yjAi 



a 2 



< 



(3.10) 



which coincides with the result of |L9| upon a trivial overall rescaling of a. The solution 
with the attractive Casimir force turns out to have a tachyonic radion, and therefore is 
unstable. 

Before concluding this subsection, we would like to briefly discuss another solution 
with non-zero U and with two positive tension branes (with a negative bulk cosmological 
constant) which was derived in |2{J. The conformal factor, for this solution, had the form 



a(t,y) = cosh 1/,2 (u;fo/) 



CO 



2 
3 



2 -"«?|A fl | 



The individual brane tensions (located at y — y\ and y = y 2 ) were given by 

Ai = — 2 tanh(o%jJ , 



(3-11) 



(3.12) 



In this case, the extremization constraint ( |2.27| ) becomes, 



25fo(A B + -f A,) 
o 



i=2 



U{byi,by 2 ) 



b=b =i 



and easily reduces to 



2A B (y 2 -yi)+ Vi 







1=1,2 



<9(fo/ 4 



-U(by 1 ,by 2 



0. 



6=6 =1 



(3.13) 



(3.14) 



As one can plainly see, the existence of this solution requires dU/db ^ 0. In [20|, this term 



was linked with Tf. 



Finally, we apply the stabilization constraint fl2.30p , which becomes 

d 2 



2u 2 (y^Ai + y 2 2 A 2 ) + ^ yi% 



i=l,2 



d(byi)d{byj 



-U 



(3.15) 



6=1 



Thus provided that d 2 U/db 2 is positive and sufficiently large, this solution is stable. 



^. Non-trivial Vi(<j>), trivial bulk potential Vb = Ab, U = 0. 

Exact solutions of this type, a massless scalar field in the bulk with non-trivial interac- 
tions on the branes (</>), were found in Refs. |22, |20[. Originally, it was hoped that such 
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a solution had some relevance for the cosmological constant problem, until more careful 
considerations [23] exemplified the necessity of fine-tuning in this model. For y\<yi< 0, 
the solution was found fOUl to be 



a (by) = smh[ujb (—y)} 



1 



2|A 



B\ 



(3.16) 



b dy sinh(ci;6 \y\) ' 

where u 2 = 8k 2 \Ab\/3, and the sign ambiguity is resolved only when the scalar potential 
on either of the two branes is specified. 

The extremization constraint ( |2.27| ) takes the form 

i=2 



Zaiyi 



5 t/2 



A B + -fV; 

6 



8 I d<p J 



0. 



(3.17) 



i=l 



and is trivially satisfied (separately for y x and y 2 ) upon the following substitutions: 

2 



V 



6|A, 



h r. 



coth 2 (c^|2/i|) , ~ ( -yj 



2\A 



B\ 



sinh (a;|yj|) 



(3.18) 



These relations allow one to express dVi/dcf) in terms of Vi and the bulk cosmological con- 
stant, which obviously requires fine-tuning: 



'dVi 



(3.19) 



Armed with these relations, we easily obtain the radion mass term from Q2.30| ), if we 
make one further assumption for d 2 Vi/d<f) 2 . In the case of linear dependence, Vi = a^ipi, the 
second line in (12.301) is zero and 



A- 



d 2 V eff 



db 2 



b=l 6 i=l,2 



lAsl-'fV 2 



(3.20) 



By using the boundary conditions (|3.18|) twice, the radion mass term may be written as 



A' 



d 2 V eff 



db 2 



-4lu\A 



B\ 



b=l 



III 



coth(o;||/i|) 2 coth(o% 2 |) 



sinh(a;||/i| 



sinh (u 1 1/2 1 



(3.21) 



The above expression is not sign definite. A simple numerical analysis reveals that the 
sign of the radion mass squared depends on the inter-brane distance (large inter-brane 
distances add a positive contribution to the radion mass term), and the distance of the 
brane system from the singularity located at yo — (small separation between the branes 
and the singularity adds a negative contribution). When \yi — 2/ 1 — \y\\ < \yi\cHti the 
stability of the system strongly depends on the location of the second brane: if it is located 
closer to the first brane rather than the singularity, the destabilizing force is minimal and 
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the solution is stable; if, however, the second brane is moved closer to the singularity, the 
radion field is inevitably tachyonic. On the other hand, for \yx\ > \yi\ cr it, the stabilizing 
force coming from the large inter- brane distance is dominant, which guarantees the overall 
stability of the system even when the second brane is placed close to the singularity. 

For non-linear interaction terms Vi, there is an extra term in the expression for the 
radion mass term proportional to d 2 Vi / dcf 2 . If the sign of this term is chosen to be positive 
then, for large inter-brane distances, this term tends to stabilize the solution even more, 
while, for small brane separations, it may reduce, or even overcome, the destabilizing force 
due to the singularity. 

As a final comment, let us stress at this point that the two-brane solution with a bulk 
scalar field, studied above, turned out to be stable, for much of the parameter space, in 
the absence of a bulk potential Vb for the scalar field. This result differs from that in [|l(| . 
where the mass of the radion is found to vanish in the limit of the vanishing bulk mass 
for the field. Our study shows that the non-trivial interaction terms of the scalar field 
with the branes are sufficient for both fixing the inter-brane distance |^U| and stabilizing 
the brane-system under small time-dependent perturbations. 



4 Conclusions 

In this paper, we have focused on the study of the stabilization of the size of extra dimen- 
sions. A variety of solutions in extra-dimensional brane-world models have appeared in the 
literature that are characterized by a constant radion field (b = 0). These solutions corre- 
spond to extrema of the corresponding radion effective potential and constraints involving 
the components of the energy-momentum tensor, for the vanishing of the first derivative of 
V e ff, have been written down for particular classes of models. 

In the context of our analysis, the derivation of those two types of constraints took place 
in the framework of the four-dimensional, effective theory for the graviton and radion field. 
We demonstrated that the variation over the size of the extra dimension reduces the integral 
constraints to a set of algebraic conditions formulated at the positions of the branes. Clearly, 
the extremization constraints give the same information as Einstein's equations. However, 
the stabilization constraint is not trivial. We were able to derive it in terms of the bulk 
and brane potentials, which may be used for a direct calculation of the radion mass term 
in many brane-world models. 

The resulting constraints were then applied for a number of static, five- dimensional 
solutions that had previously appeared in the literature. All of them were shown to satisfy 
the extremization constraint as expected, and furthermore their stability behaviour was 
studied by making use of the stabilization constraint. The two-brane Randall-Sundrum 
model was shown to be an absolute saddle point (a flat direction) of the radion potential 
while the implementation of an attractive Casimir force between the two branes led to an 
unstable configuration, in agreement with a previous work. Finally, the two-brane model 
with a bulk scalar field was shown to be stable for large inter-brane distances, i.e. when the 
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branes lie in the exponential regime of the sinh-like solution, while a potential instability 
arises for small brane separations. Of particular importance was the fact that the bulk 
potential is not a necessary ingredient for stabilizing an extra dimension with the scalar 
field. Indeed, the interaction terms between the bulk scalar field and the branes lead to a 
sufficient condition for stabilization. 
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